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1. Cosmological Phase Transitions Electroweak 

It is possible that the universe has undergone a number of phase transitions, as 
illustrated in Table 1. In most cases, it is difficult to find a signature of such a 
transition which survives to the present day. One important class of exceptions is 
when (meta)stable topological defects like cosmic stringes are formed; I will not 
deal with this important topic in the present talk. Instead I will focus on the other 
main possibility of interest, the case of a first order transition. 

A memorable example of the effect a first order phase transition could have 
was proposed by Witten [1] in 1984. He noted that when bubbles of the chiral 
symmetry broken phase form (where (q~LqR) ^ 0), baryons tend to pile up on the 
bubble walls. Neutrons diffuse quickly into the bubble interiors, but protons diffuse 
more slowly, and the spatial separation of isospin was found to have an observable 
effect on helium production in primordial nucleosynthesis. Unfortunately, lattice 
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Table 1. Some phase transitions which may have taken place during cosmo- 
logical history. 
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studies have shown that the QCD chiral phase transition is a smooth crossover 
when quarks have masses, so that in fact there is no bubble formation. 

If the QCD transition was not first order, what about the next lowest energy 
example on our list, the electroweak phase transition (EWPT)? In this case the 
bubbles would contain regions of nonvanishing Higgs field VEV, (H) ^ 0. If CP is 
violated on the bubble wall, there can be a pile-up of chiral charge, which biases 
the anomalous sphaleron interactions of the standard model to produce baryons. 
We could thus explain the baryon asymmetry of the universe. Unfortunately, the 
lattice gauge theorists have again spoiled our fun by finding that for Higgs masses 
m H > 70 — 80 GeV, the transition is a smooth crossover. (The latest limit from LEP 
is m H > 106 GeV.) However, it is much easier to change this negative conclusion 
by adding new physics (like supersymmetry) to the electroweak theory than it is 
for QCD. This will be the subject of the rest of this talk. 



2. How strong is the EWPT? 

One of the main tools for studying the strength of the EWPT analytically is the 
finite temperature effective potential, defined by 

•* i 

It is a path integral over fluctuating fields <pi, around a constant background field 
in our case the Higgs field. The fields are in imaginary time with periodic boundary 
conditions (for bosons; antiperiodic for fermions) between r = and r = [3 = 1/T. 
V c ff can be computed in perturbation theory, represented by Feynman diagrams 
like O & t one loop, and 00 or O a t two loops. The one loop term is the effect of 
a noninteracting boson or fermion gas, of which the particle masses depend on the 
background field <&, 

14-ioop =^T/^ln(l± e~ Vv^m/T^ (2) 

where the upper (lower) sign is for fermions (bosons) in the loop. This can be 
approximated in a high-temperature expansion as 

m ^ T2 (x2 for bosons) ~ $ 2 T 2 

i 

_ ^ m[g)T (bosong Qnly) ^ _ $3r 

i 

^ ! i?( 1 V +a ) +OK6/T2) (3) 

The $ 2 T 2 term is responsible for symmetry restoration at high T, while the — $> 3 T 
term gives a barrier or bump ( r - J ) in the potential, which is responsible for the 
first order transition, if it occurs. 
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Unfortunately, perturbation theory can be unreliable, especially near the phase 
transition [2]. If one starts with an arbitrarily complicated diagram contributing 
to V n -\ oop and adds an extra W boson propagator, the "cost" of the new loop is a 
multiplicative factor which parametrically has the form 



The relevant value of $ is $ c , the VEV in the broken phase at the critical tem- 
perature where V eS (0) = V eff ($ c ). If we parametrize V cS ^ AT 2 $ 2 - BTQ 3 + A$ 4 
then $ c = 2BT/X — g 3 T/X ~ (m%/m%)gT. Therefore the perturbative parameter 
e goes like X/g 2 ~ m^/mf^, and perturbation theory breaks down for heavy Higgs 
bosons, rriH > raw- 

There are several ways to combat the breakdown of pcrtubation theory at finite 
temperature. The most brute force method is to use lattice gauge simulations 
for the full 4-D theory [3]. Somewhat easier is to use dimensional reduction — 
integrating out the heavy Matsubara modes (the Fourier modes of the compactificd 
imaginary time direction, with masses m n ~ mrT) to get an effective 3-D theory 
[4,5]. The 3-D theory can then be studied on the lattice [6], much more easily than 
the 4-D theory. A third method is to compute V c g to higher order in perurbation 
theory [7-10]. This sounds unjustified, since pertubation theory was supposed to 
be breaking down, but experience shows that in fact it works rather well in the 
cases of interest — where the transition is strongly first order. 

As mentioned above, the lattice studies have established that, although there is a 
line of first order phase transitions in the T — ran plane at small tuh, it comes to an 
end (at a point where the transition is 2nd order) around wlh = 75 GeV. For larger 
ran there is no clear distinction between the unbroken and broken phases of the 
electroweak theory. For example, massive W bosons in the broken phase cannot be 
distinguished from massive composite objects (H' d^H) in the symmetric, confining 



To get a first order transition for realistic Higgs masses, we need to add new 
physics which couples significantly to the Higgs boson. Let's see how supersymme- 
try can do this. 

3. Adding Supersymmetry 

Recall that at one loop, it is the cubic term, —Tra 3 (H)/\2n that gives bump in 
the potential (r\J ) hence a first order transition. In the minimal supcrsymmetric 
standard model (MSSM), we have two Higgs doublets, and the top squarks (?l and 
is) couple strongly to the second one, H 2 - Ignoring small terms involving the weak 
gauge coupling g, the t L -t R mass matrix has the form 



We need at least one of the SUSY breaking soft masses, ra\j or TOq, to be small 
so that there will be an eigenvalue ra{H) which is really cubic, <~ H 3 , not (m 2 , + 





phase. 
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y 2 H 2 ) 3 / 2 , since the latter form does not give a true bump in the potential [8]. On the 
other hand, one of mfj or m,Q should be large so that the stop radiative correction 
to the Higgs mass can be big enough to satisfy the experimental constraint: 



O 



mf , I m; m; 



^ln 



9 I 2 

v V mi 



> (96 GeV) 2 (6) 



(The limit on m H in the MSSM is about 10 GeV weaker than in the SM.) The 
precision electroweak p parameter (a.k.a. ei) dictates that the iL-like squark should 
be the heavy one, hence the f^-like squark is light. 

It turns out that two-loop effects are crucial for getting a strong enough phase 
transition[7-10]. Diagrams like ©-, with a gluon or Higgs boson as one of the 
internal lines, contribute a term 

AV ff = ~(8g 2 s - 3y 2 sm 2 p)T 2 m 2 (H)ln('"'>" ! ' 



T 

-CT 2 H 2 \n(^\ (7) 



to I4ff , whose form is unlike any generated at one loop. One can show that such a 
term shifts the critical value of the Higgs VEV according to 




(8) 

where B is the part coming from the cubic term. (H)/T is the relevant measure of 
the strength of the transition, as we will discuss below. 

In ref. [10] we have performed a Monte Carlo search of the MSSM parameter 
space to find those which give a strong enough phase transition for electroweak 
baryogenesis. We can summarize the resulting constraints on the squark masses as 
follows: 

120 GeV < m~ tR < m top (9) 
m iL > 265 GeV x e (™«- 95 GcV)/9.2 ^ 

It is often said that electroweak baryogenesis puts an upper limit on the Higgs 
boson mass, but this is not correct in the MSSM, where the light right stop is doing 
most of the job of making the transition first order. However, as we can see from 
eq. (ft0|), it is true that the left stop quickly becomes unnaturally heavy as m# 
is increased. Thus we are pushed to a rather strange corner of parameter space, 
where tn is extremely light, and ti extremely heavy. If mu turns out to be much 
heavier than its current experimental limit, there must be some other new physics 
accounting for its mass. 

More promising is the prediction (^) for the stop mass. Although Tevatron 
searches for the decay t — > Xi c are on ly beginning to probe the region of interest 
(fig. 0(a)) [11], in Run II this decay will probe up to the top mass if m^o is in the 

right range (dark region of fig. [l](b)), and the decays t — > bxt or t — > bWxi ma y 
be more revealing, if the chargino is light (lighter regions) [12]. 
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CDF Sensitivity to stop in Run II 




M(t,) (GeV/c ) Mstop (Gev/c") 



Figure 1. (a) Present Tevatron limits [11] and (b) sensitivity in Run II for 
the stop mass [12] in the (neutralino mass)-(stop mass) plane. 

An interesting consequence of such a light stop is that its bare (mass) 2 must be 
negative. In the absence of left-right mixing, m~ = TO^ + m 2 , so < m t implies 

to 2 / < 0. This can cause an instability toward condensation of the tR field in the 
early universe, when (H) is still zero, which would break SU(3) co i or [9]. Indeed, 
the lattice study of ref. [13] finds a phase diagram similar to fig. ||(a), which shows 
that color-breaking occurs when (— m 2 ;) 1 / 2 exceeds 60 — 70 GeV, depending on 
the critical temperature T c . In ref. [14] we have constructed the two loop effective 
potential V e s(tn,H) for stop and Higgs fields, and studied the possibility that the 
universe might temporarily enter the color-breaking phase, before finally tunneling 
to the EW-breaking vacuum which we inhabit now. We find that because of the 
potential barrier separating the two minima, shown in fig. |^(b), the rate of tunneling 
is so small that if the universe ever enters the color-breaking minimum, it stays there 
forever. This conclusion could however change in the presence of R-parity violating 
interactions like y332At R b R s c R e a b c which would lower the barrier. 

4. Baryon Asymmetry of the Universe and Baryogenesis 

I will now review some basics about baryogenesis and describe some recent develop- 
ments, before giving the latest details on electroweak baryogenesis in the MSSM. It 
is very unlikely that we live in a baryon-antibaryon symmetric universe since there 
is no evidence of antigalaxies colliding with galaxies, which would give an intense 
source of gamma rays. If there are regions consisting of antimatter outside of our 
Hubble volume, it is difficult to imagine how separation from ordinary matter could 
occur on such a large distance scale. Big Bang nucleosynthesis tells us that 
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Figure 2. (a) Phase diagram and (b) effective potential, V e g {Ir, H), for color 
breaking due to stop condensation. 

2 < ^o^^-^xlO- 10 < 3 (11) 

n 7 

Yet in the early universe the most natural initial condition is equal numbers of 
baryons and antibaryons, since there are many possible B-violating interactions 
that could be in thermal equilibrium at high temperatures, for example a dimension 
9 operator like (udd) 2 /A 5 , which would cause neutron- antincutron oscillations, or 
the interactions with heavy X gauge bosons in grand unified theories. In fact 
we need not look so far afield, since sphalerons violate baryon number within the 
standard model itself! 

It is well known that Sakharov's three conditions must be fulfilled to generate 
a baryon asymmetry: (1) baryon number violation; (2) C and CP violation; (3) 
loss of thermal equilibrium for the B-violating interactions. The first and second 
conditions are present within the SM, but (2) is too weak for baryogenesis, and (3) 
is not fulfilled at all. We have already seen how the MSSM can increase (H) /T; this 
is what is needed to make the sphalerons go out of equilibrium inside the bubbles 
that form during the EWPT. The MSSM can also cure the problem of getting 
strong enough CP violation, as we will discuss in the next section. Here, I will only 
mention some of the other proposals for baryogenesis which are new or currently 
popular. 

Baryogenesis via leptogenesis [15] is one of the most plausible alternatives to 
electroweak baryogenesis. In analogy to GUT baryogenesis, heavy sterile neutrinos 
decay out of equilibrium, producing a lepton asymmetry, which is converted by 
sphalerons into the baryon asymmetry. The predictions can be related to neutrino 
masses in a GUT framework like SO(10) [16]. 

The Affleck-Dine mechanism [17] has long been one of the most efficient baryo- 
genesis mechanisms. It uses the fact that SUSY scalar potentials from D-terms can 
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often have flat directions which generate a huge baryon number when the flatness 
is lifted and the field evolves by spiraling in the complex plane. It was recently 
pointed out that this can be combined with leptogenesis to give a minimal super- 
symmetric baryogenesis model in which the flat direction is a linear combination 
of H 2 (the second Higgs doublet) and L e (the selectron), using physics which is 
already needed for generating neutrino masses [18]. 

Large extra dimensions can present a serious challenge to baryogenesis by con- 
straining the reheat temperature after inflation to be very low [19]. The Randall- 
Sundrum alternative of warped compactification [20] evades this problem; so per- 
haps do intrinsically "braney" approaches to baryogenesis [21]. 

Other novel ideas make use of the phase transition in left-right symmetric models 
[22] and decaying primordial black holes [23]. 

5. Electroweak Baryogenesis in the MSSM 

Electroweak Baryogenesis in the MSSM is one of the most carefully studied ideas for 
baryogenesis, owing to its close ties to present-day phenomenology and accelerator 
searches [24-29]. Its basic mechanism [30] is intuitively easy to understand: par- 
ticles interact in a CP-violating manner with bubble walls, which form during the 
first order electroweak phase transition, when the temperature of the universe was 
near T = 100 GeV. This causes a buildup of a left-handed quark density in excess 
of that of the corresponding antiquarks, and an equal and opposite right-handed 
asymmetry, so that there is initially no net baryon number. The left-handed quark 
asymmetry biases anomalous sphaleron interactions, present within the standard 
model, to violate baryon number preferentially to create a net quark density. The 
resulting baryon asymmetry of the universe (BAU) soon falls inside the interiors 
of the expanding bubbles, where the sphaleron interactions are shut off (provided 
that (H) > T), and thus baryon number is safe from subsequent sphaleron- induced 
relaxation to zero. 

However, the correct way to treat the generation of the chiral quark asymmetry 
in front of the bubble wall is controversial. In the simplest model, the top quark 
has a spatially varying complex mass, m(x) = \m{x)\e iet ^ x \ which gives rise to CP- 
violating quantum mechanical reflection of quarks as they pass through the bubble 
wall. It also induces a CP-violating classical force on the quarks. In rcf. [31] it 
was shown that the latter is the more appropriate treatment when the bubble wall 
thickness l w is large compared to the inverse temperature, which is the case in 
the MSSM: l w = (10 ± 4)/T [29]. Furthermore it is possible to rigorously derive 
the way in which the CP violating force influences the particle transport, using 
the Boltzmann equation; no such complete derivation yet exists in the quantum 
reflection formalism. 

In the MSSM, the top quark mass does not have a CP-violating phase; however 
the charginos do; their mass matrix has the form 

faM* _ K, S+) R ( ™ {x) **<«> ) ( 1 1 ^ + h.c. (1 2 > 
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Since the SUSY parameters /i and m,2 can be complex, the mass eigenstates can 
have spatially varying complex phases in the wall. There thus arises a classical 
force which separates the two kinds of Higgsinos, h\ L and h\ L , in front of the 
wall. This kind of asymmetry is not enough to bias the sphaleron interactions, but 
scatterings, such as hig — > £l£r, will partially convert the Higgsino asymmetry into 
a chiral quark asymmetry, n qL . Once the latter is determined, it is straightforward 
to integrate the rate of baryon violation by sphalerons, governed by the equation 

(13) 

The rate of sphaleron interactions per unit volume has been measured by lattice 
simulations to be T sph = (20 ± 2)a 5 w T 4 [32]. 

To determine n qL , we must solve a set of coupled diffusion equations for the 
various species i of particles in the plasma. They have the form 

- Ai" - v w rii + T ijk (m - nj - n k ) = Si, (14) 

where Di is a diffusion coefficient (of order the inverse mean free path), v w is the 
bubble wall velocity, Tijk is the rate of interactions of the type i — > j + k (given 
as an example), and Si is a CP-violating source term arising from the force on the 
particles or from quantum reflections. It is possible to show that the source term 
in the Higgsino diffusion equation is related to the force F by the thermal averages 

S(x) = -^-{v x F(x))', (15) 

and, by solving the Dirac equation in the WKB approximation, the CP-violating 
part of the force is related to the complex Higgsino mass me by F — 
(s/2E 2 )(m 2 6'Y , where s — ±1 is the spin. Moreover the combination m 2 9' is 
given by 



m 2 £ 



flHm^) {HiH , +H[H2h (16) 



2(m3_ — m_) 



where m| are the two eigenstates of the mass matrix in (12). A remarkable feature 
of this expression is the relative + sign between HiH' 2 and H' l H2- Other authors 
[24,26] had previously obtained a source proportional to H\H' 2 minus H[H2, which 
is highly suppressed, because H1/H2 tends to be constant within the wall, to within 
a part in 10 2 or 10 3 [33,10]. The origin of the discrepancy is that (for technical 
reasons) the previous authors considered only the linear combination of Higgsino 
densities n% — ni , whereas our classical force is providing a source for + n~ h 

By solving the difffusion equations and numerically integrating the baryon viola- 
tion rate equation, we obtained the baryon asymmetry as a function of the model- 
dependent parameters of the MSSM [29] . Figure ||(a) shows how 7710 varies with 
the bubble wall velocity v w in a typical case, for the allowed range of values of the 
bubble wall thickness, where we assumed maximal CP violation, Im(m2yu) = \m2n\- 
Since typically 7710 is ~ 1000, but we only need (2 — 3), we see that the CP-violating 
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phase need not be maximal, but could be (2 — 3) x 10~ 3 . This is good news, since 
the neutron and electric dipole moment searches give constraints which are of this 
order, unless some kind of fine tuning is invoked. Baryon production peaks at small 
wall velocities near v w = 1CP 2 . Interestingly, recent estimates of v w in the MSSM 
give values which are this small [34]. Figure ||(b) shows the contours of constant 
CP phase which yield ry 10 = 3 in the chargino mass parameter plane. 



6. Conclusions 



In this talk I have discussed one of the main "applications" of cosmological phase 
transitions, baryogenesis. Although there are many imaginative ideas for getting 
the baryon asymmetry, the most mainstream ones are those that are most closely 
related to phenomenology. Baryogenesis from leptogenesis is appealing because it 
can potentially make contact with neutrino masses. Electroweak baryogenesis is in- 
directly testable by searches for Higgs bosons, the top squark, charginos, neutralinos 
and electric dipole moments. We have seen that it is relatively easy to generate a 
large enough baryon asymmetry in the electroweak model with the MSSM, but it is 
not so easy to get a strong enough phase transition to safeguard it against washout 
by sphalerons: a large hierarchy between the left and right stop masses is needed. 
It would be nice to have some more robust way of strengthening the transition. 
One possibility is adding a singlet Higgs field [35]. An interesting proposal is to re- 
heat after inflation to temperatures below the electroweak transition, but rely upon 
nonequilibrium effects, similar to parametric resonance, to produce sphalerons [36]. 
Another idea is to modify the expansion rate of the universe instead. It has been 
noted that extra dimensions could have this effect at sufficiently early times [37]. 
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